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Abstract 

The Cabibbo-Kobayashi-Maskawa (CKM) parameters are investigated in detail in recent pre- 
dictive models which are based on low-energy non-abelian discrete family symmetries. Some of 
the models can already be excluded at the present precision of the determination of the CKM 
parameters, while some of them seem to survive. We find that to make the uncertainties of the 
theoretical values comparable with the assumed uncertainties of ~ 1° and ~ 2° in (j)2{oi) and 03(7), 
respectively, at about 50 inverse atto barn achieved at a future B factory, it is necessary to reduce 
the uncertainties in the quark masses, especially that of the strange quark mass by more than 60%. 

PACS numbers: 12.60.Jv,11.30.Hv, 12.15.Ff, 14.60.Pq, 02.20.Df 
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I. INTRODUCTION 



The success of the standard model (SM) suggests that we are very close to a more 
fundamental theory for elementary particle physics. Yet, we do not know how the SM 
should be extended, except that the mass of neutrinos with their mixing and a dark matter 
candidate have to be incorporated in the extension. The Higgs sector of the SM indicates that 
the extension may take place around TeV scale, and supersymmetry is widely believed to be 
the best candidate to increase the natural energy scale of the theory. Yet another problem 
is the Yukawa sector, because the most of the free parameters of the SM are involved there 
and the SM does not provide with a principle how to fix its structure. Moreover, simple 
supersymmetry does not soften the problem of the Yukawa sector. 

A natural way to provide with a principle for the Yukawa sector is the introduction of a 
family symmetry. A family symmetry is not necessarily adequate to explain the observed 
hierarchy of the fermion masses. It can however relate the fermion masses and mixing 
parameters. That is, mixing parameters may be related to mass ratios. Note that the 
classic relations such as sin^c* — V^^v/"^ iKtb/Ktl — \fnh^rnc ^, E| jsl had 

not been derived from a family symmetry. Recently there have been a growing number of 
interests in family symmetries. Most of the recent papers deal with the large neutrino mixing 
(see for instance [30]), because a large mixing may be associated with a family symmetry. 
As for the quark masses and mixing, tremendous works on their ansatz have been done (see 



for instance [lOj). However, there is an almost no-go theorem [ll| saying that there exists 
no viable low-energy family symmetry in the SM to understand the fermion mass matrices. 
Therefore, if the fermion mass matrices should be derived from a family symmetry, one has 
to extend the SM. One of the possibilities is to extend the Higgs sector such that it also 
forms a family |l2[| . In fact a renewed interest in this approach to the quark sector has been 



recently aroused |l3|-|55| 



In this paper we are interested in predictive flavor models with a low-energy discrete 



family symmetry that are testable at future B factories such as SuperKEKB 56| and Super 



Flavour Factory [57j. We met a set of the following selection criteria for the models: 

1) The family symmetry should be a low-energy symmetry. That is, we do not include here 
family symmetries at GUT scale^. 

2) The family symmetry should not be hardly broken. If it is hardly broken, there is in 
general no quantitative prediction of the symmetry. So, we do not consider textures which 
are not supported by a symmetry. 

3) The model should describe 10 observables, six quark masses and four Cabibbo-Kobayashi- 
Maskawa (CKM) parameters, by less than 10 parameters. 

4) The model should be renormalizable. 



To our knowledge there are only four models proposed in 20l-l22| that satisfy the selection 



^ Recent GUT models with a family symmetry can be found in [43|-|55j. 
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criteria above. 

In the next section we will start by summarizing the present values of the CKM pa- 
rameters and quark masses. We will find that the uncertainties in the light quark masses 
have to be largely reduced to make the uncertainties of the theoretical values of the CKM 
parameters comparable with the assumed experimental uncertainties of future B factories. 



The lattice calculation [58| is indeed an on-going project to reduce the uncertainties in the 
igh t quark masses. In this paper, however, we will not use the quark mass values given in 
58|, because the uncertainties due to the absence of the sea strange quarks have not been 

included. 



II. QUARK MASSES AND CKM PARAMETERS 

Our strategy to compare the theoretical values with the experimental ones is to use the 
six quark masses along with two of the CKM parameters to plot the theoretical values in 
the plane defined by other two CKM parameters. From the quark masses in the MS scheme 
given in Particle Data Group 2009 |59| we obtain the quark masses at /i = Mz 



60 



mu{Mz) = 0.98 ~ 2.15 MeV , md(Mz) = 2.28 ~ 3.90 MeV, 

ms{Mz) = 45.5 ~ 84.5 MeV , m^iMz) = 0.65 ~ 0.75 GeV, (1) 

mb{Mz) = 2.85 ~ 3.02 GeV , mt(M^) = 179.2 ~ 183.7 GeV. 

The input parameters we shall use are the ratios at /i = Mz, i.e. 

mjmt = (0.54 ~ 1.18) x 10"^ , mjrrit = (0.36 ~ 0.41) x 10-^ 

md/rrib = (0.79 ~ 1.34) x 10^^ , m./mf, = (0.16 ~ 0.29) x 10"\ (2) 

Given the ratios above, there are two independent mass scales, one for the up quarks 
and the other for the down quarks: if necessary we use the fixed values mb{Mz) = 



2.9 GeV , mt{Mz) = 181.4 GeV. We further require that |59| 



ms/nid = 17 ~ 22 , !!^!i±Z!M = 2.5 ~ 5.5 MeV , ^""^ = 25 ~ 30 , 

2 mu + nid 

mu/rrid = 0.35 ~ 0.6 , — — ^ = 30 - 50 (3) 

rrid - ruu 

are satisfied. We also use 

\Vus\ = 0.2236 ~ 0.2274 , \Vcb\ = 0.0401 ~ 0.0423 (4) 

as the two input parameters of the CKM parameters. Then the theoretical values are 
compared with two sets of the experimental values: 

PDG:i59| 
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TABLE I: The S3 and Z2 assignment in the model I. 2ll | 



and 



p = 0.1425 ± 0.0235 , f] = 0.348 ± 0.016, 

02(a) = (88.5 ± 5.5)° , ^3(7) = (76 ± 31)° , sin20i(/3) = (0.681 ± 0.025), 
\Vub\ = (3.93 ± 0.36) X 10"^ , \Vtd/Vts\ = 0.209 ± 0.007, 



(5) 



UTFit:|61| 



p = 0.154 ± 0.022 , f] = 0.342 ± 0.014, 

02(«) = (92.0 ± 3.4)° , 03(7) = (65.6 ± 3.3)° , sin20i(/3) = 0.695 ± 0.020, (6) 
\Vub\ = (3.60 ± 0.12) X 10"^ , \Vtd/Vts\ = 0.209 ± 0.0075. 

III. THE MODELS AND THEIR CKM PARAMETERS 



In this section we consider four different flavor models with three different family sym- 
metries. Except the model II of j22| these models should be supersymmetric to obtain 
desired mass matrices for the quarks. The family symmetry of the models III and IV of 



20| extends to the leptonic sector so that there are nontrivial, testable predictions in that 
sector, too^. The SU{2)l doublets of the quarks and Higgs bosons are denoted by Q and 
}ju,d^ respectively. Similarly, SU{2)l singlets of the quarks are denoted by and D^. In 
the supersymmetric models they are superfields, and they are ordinary fields in the non- 
supersymmetric case. 



A. The model I |2l|] 



The first model, which is proposed in 21|, is the supersymmetric model with an 5*3 x Z2 
family symmetry. The S3 x Z2 assignment is given in Table HI and the S3 x Z2 invariant 



See 



3,0 for an alternative assignment of the leptons to obtain the maximal mixing of the atmospheric 



neutrinos. 
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cubic superpotential for the Yukawa interactions in the quark sector is given by 
W, = Y^QsU^H!^ + Y,-{Q,H^ + Q^H^m 

+Y^Q,DIH', + Y^\Q,Dl + Q^DDH^ + Yf [Q.DlHi + Q^DlMf) , (7) 

where the subscripts 1, 2 and 3 stand for the two components of the 5*3 doublet and for the 
5*3 singlet, respectively. After the spontaneous electroweak symmetry breaking, the Higgs 



I ^0. 




FIG. 1: The theoretical values of the model I at 90% CL in the p—rj plane. Two experimental values ([5]) 
and ([6]) are also plotted. The best fit point is denoted as x . 



doublets should acquire vacuum expectation values (VEVs) < H^''^ >= v^''^. Then the 
quark mass matrices are'^ 



(8) 
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Note that 



(9) 



is assumed, and supersymmetry makes it possible to obtain this 



relation naturally 21 1. 

All the elements of these matrices can be made real and positive by an appropriate 
redefinition of the quark fields. Therefore, the real quark mass matrices can be diagonalized 
by orthogonal matrices as follows. 



tj 5 



O^MuM^Ou = disig{ml,ml,m 
O^MaMjOd = diagK, m2,m2) , 



(10) 
(11) 



This model is very similar to the ansatz considered in [62| . 
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TABLE II: Dj assignments of the matter fields, 
and the CKM matrix is written as 

V = OlPO, , 



(12) 



where P = diag(l, e*''^, 1). There are only eight independent parameters, so that we can 
calculate two remaining physical quantities. 

The CKM matrix V can be approximately obtained in a closed form, and one finds 
for instance, 



21 



ub 



Vr 



cb 



XW 



m. 



u 



td 



~ — ~ W , — ^ -tanod ~ — . 

yz V "^c Vts V m.^ 



rrid 



(13) 



Using the values given in ([T]), ([2]) and (jl]), we find 

(1.4 ~ 2.4) X 10~^ 



Vub\ 

Vtd 



IK 



cb\ 



nir. 



ts 



rrid 




0.213 ~ 0.242 



(14) 
(15) 



As we can see from ([5]) and ([6]), the ratio |Vfd/Vis| is consistent with the experimental value, 
while \Vub\ is not. We have performed a systematic, numerical analysis of the theoretical 
values in various planes. Fig. [1] shows the case of the model I at 90% CL in the p—f] plane. 
We clearly see that the model I is not consistent with the experimental observations^. 



B. The model II |22|] 



Next we consider the Dj flavor symmetric model of 22|, where a supersymmetric exten- 
sion is not always necessary in this model. The D-i assignment is given in Table [TTl and the 
Dy invariant Yukawa Lagrangian of the quark sector can be written as 



+Y^{Qid2 + Q2di)Hl + Y,\QiH2 + Qai^i )4 
+Y^^Qs{Hid2 + H^d^) + YfQ.H^ds + h.c. , 



(16) 



^ Changing the Z2 assignment appropriately, one can interchange the mass matrices for the up and down 
quarks. In this case, however, one obtains a negative p [2l|, which is excluded experimentally. 
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which yields the quark mass matrices of the form 



\ 

Ffe" J 



f K" 
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/ 



yd Q yd „,d 



d „,d yd 
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d „,d I 
'^3 / 





V 
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(17) 



(18) 



where < H^''^ >= v^''^. Except for ^, the parameters a, b, c and d of can be made real and 
positive. That is, there are nine independent parameters in the quark sector, and we can 
calculate one physical quantity. The mass matrix of the down quarks can be diagonalized 




FIG. 2: The theoretical values of the model II at 90% CL in the p— 77 plane. Different colors indicate 
different intervals of the mass ratio ms/md'- 17 ^ 18 (red), 18 ^ 19 (green), 19 ~ 20 (blue), 20 ~ 21 
(orange), and 21 ~ 22 (brown). Two experimental values ^ and ^ are also plotted. The best fit point of 
the model II is denoted as x . 



by a unitary matrix Vd as 

and the CKM matrix is then given by 

V = Vd, 



(19) 



(20) 



because the mass matrix of the up quarks Mu is diagonal, which is a consequence of the 
family symmetry. If one assumes that -C b^ and ^ 1, one finds various approximate 
relations such as 



221: 



nib — y/c^ + d'^ , nis — 
bd 
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(21) 
(22) 
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TABLE III: The Qq assignment of the matter multiplets. The assignment is the same for the model 
III and IV. 



Instead of using these approximate analytic expressions, we have performed systematic nu- 
merical analyses. The theoretical values in the p—f] plane are plotted in Fig. O Different 
colors mean different intervals of the mass ratio ms/m^: 17 ~ 18 (red), 18 ~ 19 (green), 
19 ~ 20 (blue), 20 ~ 21 (orange), and 21 ~ 22 (brown). The best fit for the set of parame- 
ters, sin20i(/3),03(7),mrf/mfe,m,/m{„ \Vcb\ is given by 

sin2(/)i(/3) = 0.682, ^3(7) = 62.314°, |K.| = 0.2256, l^bl = 0.0411, 

ma/mb = 1.120 x 10"^ m,/mfe = 0.191 x 10"^ (23) 

with x^/dof = 0.51. As we can see from the figure, the Dj model is consistent with the 
experimental observations. More precise measurements of the CKM parameters as well as 
more precise determinations of ms/m^ are needed to confirm or exclude the model. If the 
mass ratio m^/mrf turns out to be larger than 18, the model may run into problems. 



IV. THE MODEL III [20|] AND IV 



The last two supersymmetric models are based on a Qq family symmetry 20|]. In contrast 
to the previous two models, the family symmetry of these models extends to the lepton 
sector. The Qq assignment of the leptons given in j25| indeed leads to the maximal mixing 
of the atmospheric neutrinos. Note, however, the parameter space of the model IV has not 
been discussed previously. The finite group Qq allows complex representations, and the Qq 
assignment of the matter multiplets is given in Table lllli 

The superpotential for the Yukawa interactions in the quark sector is given by 

+Y:{Q,U:i + Q,Ut)Hl 

+Y^Q,HtDl + Y,\Q,Hi + Q2H'i)Dl + Y^^Q,{HtDl - H^Dl) 
+Y^{Q^D^^ + Q2Dl)Hi . (24) 

To make the model predictive there are two crucial requirements: (1) the VEV alignment 

. ttuA ^ , TT^A ^ u,d , rru.d ^ u,d 

which can be achieved by an accidental permutation symmetry H^''^ ^ iJg i'^ the Higgs 
sector, and (2) CP is spontaneously broken. The second requirement can be relaxed to that 
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the Yukawa couplings are real without contradicting renormalizability^. Then the quark 
mass matrices can be written as 







'^3 ^1 

y^" 

yu yu yu ^ 



yu 
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„,d 



yd yd ^d \ 



'd „,d vd 

Y^ vf 

'd „,d vd „,d vd „,d 



\ -Y^ vf Y^ vf Y^ vl j 



(26) 



with complex VEVs. By making an overall 45° rotation of the doublets Q, and D"^ 



PDG 




0.1 0.1 _ 0.3 

P 



FIG. 3: The theoretical values of the model III at 90% CL in the p—r\ plane . Two experimental values ([5| 
and ([6]) are also plotted. The best fit point is denoted as x . 

in the space of the family group, we obtain nearest neighbor interaction (NNI) type mass 
matrices: 



/ c„ \ 

-Cu hu 
\ du tu) 



( Crf \ 

-Cd hd 
\ dd Cd ) 



(27) 



All the elements of these matrices can be made real by a suitable redefinition of the quark 
fields. Then the real matrices can be diagonalized by orthogonal matrices as 



OlMuM^Ou = diag(m2, m2,m2), 
O^MdMjOd = dia.g{ml ml, ml), 

and the CKM matrix takes the form 

V = OlPOd. 



(28) 
(29) 



(30) 



where P = diag(l, e 



2ie Ae\ 



The phase rotation matrix P has only one angle 9, which is the 



consequence of a spontaneously and softly broken CP. 



40( 1 that to trigger complex VEVs with the minimal content of the chiral supermultiplcts 



^ It has been found 

given in Table IIIIl the family symmetry and CP should be at least softly broken. The most economic 
breaking can be achieved by the b-terms in the soft-supersymmetry breaking sector. 
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The full nine dimensional space of the parameters can be divided into two non-equivalent 
regions. The difference between the two regions, the region III and IV, can be reduced to the 
sign of Q, so that without loss of generality the other parameters, i.e. bu,Cu. ■ ■ , in (j^Tj), 
can be assumed to be positive real numbers. (The parameter space for the model IV has 
not been considered in 20|.) Accordingly, we define two models; the model III for positive 
Cd and IV for negative c^. 




FIG. 4: The theoretical values of the model IV at 90% CL in the p-i] plane (left) and \Vtd/Vts\~(j)3{j) 
plane (right). Two experimental values ^ and ([5]) are also plotted. The best fit point is denoted as x. 




FIG. 5: The theoretical values of the model IV at 90% CL in the 03(7)— 02 (a) plane, where two predicted 
regions for ?«,/mb = 0.0220±0.0007 and md/mb = (1.200±0.028)10-3 (red) and for mjmb ^ 0.0170±0.0007 
and md/nib = (0.930 ± 0.028)10"'^ (blue) are plotted. Two experimental values ([5]) and (|6]) are also plotted. 
The best fit point is denoted as x . 



According to 63|, the CKM matrix can be approximately written in a closed form. One 



10 



finds, for instance, 



Vus ^ -VdJ — + yuJ — e''" , 31 
^' nis v rrir 



y ^ Vd ^ g2i0 _ ^« !^ giff C32) 



Vub ^ 




yd V 



( - . 1 , (33) 



wfiere j/^j = l/c„ and = For ~ 1 we can reproduce tlie classic relations of [l|-|3| 

IK.I - ./^ + 0(m„/m,). (34) 

We have performed numerical analyses on the CKM parameters of the models in detail, 
some of which are presented in Figs. [3]-l5l where Fig. [3] is that of the model III, and Figs. H] 
and [5] are those of the model IV. 

The best fit for the model III is given by 

sin20i(/3) = 0.664, 02(«) = 81.397°, |Ks| = 0.2250, IK^I = 0.0416, 

m^jmt = 0.860 x 10~^ mjmt = 0.375 x 10"^ (35) 
TTid/mh = 1.065 X 10"^ ms/rrih = 0.199 x 10^^ 

with = 2.68 for the set of parameters sin20i(/3), (l)2{a),mu/mt, rric/mt, md/rrih, nis/mi,, \ Vus\, \ Vcb\ 
As we can see from Fig. [3l the model could be excluded; the PDG value in the p — f] 



plane 59| is about 2a away from the theoretical values. A more precise experimental 



determination of 02 (a) will be more supporting the conclusion. 



The model IV, on the other hand, seems to be consistent with the PDG values 59 1 



as well as with the Utfit group values [6l|. In Fig. \5\we plot two predicted regions for 
mjmb = 0.0220 ± 0.0007 , ma/mb = (1.200 ± 0.028)10-3 (red) and m,/mfe = 0.0170 ± 
0.0007 , md/nib = (0.930 ± 0.028) 10^^ (blue), which should be compared with We see 
that precise measurements of 02(«) can distinguish the two regions. We also find that the 
smaller the is, the larger is 02 («)^- 

The best fit for the set of parameters of the model IV, 



^ While a smaher ms{Mz) has been recently reported in [6J|, the lattice calculation of [58[ including elec- 
tromagnetic interactions and non-perturbative renormalization effects leads to a slightly larger ms(2GeV) 
than that of 01. 
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p,T],mu/mt,mc/mt,md/mb,ms/mb, |Ks|, is given by: 

p = 0.144, f] = 0.352, \Vus\ = 0.2254, = 0.0411, 
niu/mt = 0.704 x 10"^ mc/mt = 0.390 x 10"^ (36) 
ma/mb = 0.965 x 10"^ m,/mb = 0.177 x 10"^ 

with = 1-03, and 

(/,3(7) = 66.698°, (/)2(a) = 87.650°, |K.| = 0.2254, iKbl = 0.0411, 
m^/mt = 0.788 x 10^^ mjmt = 0.390 x 10^^ (37) 
ma/mb = 1.075 x 10"^ m,/m6 = 0.194 x 10"^ 

with = 0.44 for the set of parameters 03(7), (j)2{o:),mu/mt, mc/rrit, md/rrib, nis/mb, \ Vus\, \ Vcb\ 
V. CONCLUSION 

We have investigated four predictive flavor models that are based on a low-energy discrete 
family symmetry. Two of them may be excluded by the present precision of the measure- 
ments of the CKM parameters. The problem of the model I has been already pointed out 



m 



2l|. For the model II a more precise determination of the quark masses is very crucial, 
especially the mass ratio nig /ma. We have found that the smaller the mass ratio is, the 
better is the chance for the model II. We have also found that to confirm or exclude the 
model IV, more precise determinations of 03(7) and 02(tt) as well as of the quark masses 
are indispensable. The uncertainty in the strange quark mass should be less than a little 
more than 10% to make it comparable with the assumed uncertainties of ~ 1° and ~ 2° in 
02(«) and 03(7), respectively, at about 50 inverse atto barn achieved at a future B factory. 
In this letter we have restricted ourselves to low-energy predictive models. There are also 



high-energy predictive models [43|]- 55|, which will be the target of our future investigations. 
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